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ABSTRACT 

We describe recent Chandra ACIS observations of the Vela-hke pulsar PSR 
J2021+3651 and its pulsar wind nebula (PWN). This 'Dragonfly Nebula' displays 
an axisymmetric morphology, with bright inner jets, a double-ridged inner nebula, 
and a ~ 30" polar jet. The PWN is embedded in faint diffuse emission: a 
bow shock-like structure with standoff ~ 1' brackets the pulsar to the east and 
emission trails off westward for 3-4'. Thermal (kToo = 0.16 ±0.02 keV) and power 
law emission are detected from the pulsar. The nebular X-rays show spectral 
steepening from P = 1.5 in the equatorial torus to P = 1.9 in the outer nebula, 
suggesting synchrotron burn-off. A fit to the 'Dragonfly' structure suggests a 
large (86±1°) inclination with a double equatorial torus. Vela is currently the 
only other PWN showing such double structure. The >12kpc distance implied 
by the pulsar dispersion measure is not supported by the X-ray data; spectral, 
scale and efficiency arguments suggest a more modest 3-4 kpc. 

Subject headings: pulsars: individual (PSR J2021+3651) — stars: neutron — 
gamma-rays — observations 

1. Introduction 

Pulsars inject a relativistic particle wind into their surroundings, tapping their rotational 
kinetic energy. This energetic wind is concentrated toward the spin equator and, shocking 
against the surrounding medium, causes particle acceleration and pitch angle sca ttering that 



results in a synchrotron nebula (IRees fc Gunru Il974l : iKennel fc Coronitil Il984l ). Approxi 



mately 50 pulsars have an associ ated PWN, and the polar j et /equatorial torus morphology 



is now seen to be quite common (iKargaltsev &: Pavlov! l2008l ). Many of the mechanisms gov- 



erning the energetics and evolution of pulsar/PWN systems are still unknown, such as how 
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rotational kinetic energy is converted to particle outflow, the manner in which polar jets 
are confined, and the details of charge acceleration. The structure of PWNe offer insights 
into these puzzles, and also probe pulsar geometry, relativistic shocks, and the properties 
of the ambient medium. PSR J2021+3651 (hereafter J2021) is interesting for its highly 
structured PWN, and is also a likely counterpart of the EGRET gamma ray source GEV 
2020+3658. Pulsars detected in GeV gamma rays are currently of great interest with the 
successful launch of AGILE and the imminent launch of GLAST. 

Timing observations of this 103.7 ms pulsar show that J2021 is young and energetic 
[characteristic age Tc = P/{2P) = 1.7 x 10^ yr and spindown luminosity E = An'^IPP'^ = 
3.4 X 10 ^^ ergs s~^ ], despite being rather radio faint with a 1.4 GHz fiux density Suoo ~ 
0.1 mJy jRoberts et aPbooi iHessels et al.l 120041 ). The dispersion measure of 369 cm ^ pc 
places it ~12 kpc away on th e far edge of the oute r spiral arm, according to the NE2001 
electron model of the galaxy (jCordes fc Lazidl2002l ). Alth ough its spin charac teristics are 
similar to those of other 7-ray pulsars [Vela: d ^ 300 pc (IDodson et al.ll2003l ). P = 89.3 



ms, E = 6.9x10^^ ergss"^ B1706-44: d ~3 kpc (IRomani et al. 



20051). P = 103 ms, E 



3.4x10^^ ergss ^], the 12 kpc distance assume d from the dispersion measure would make it 
extraordinarily 7-ray efficient, as discussed by [Roberts et al.l (120021 ). 



Roberts. Romani. &: Kawai I (I2OOII ) noted the possible connection between GEV J2020+3658 



and ASCA source AX J2021. 1+3651. The ASCA data revealed spatially unresolved emission 
extending ~8' east to west and ~4' north to south, with a hydrogen column density of A^h 
= 5.0+2.5 X lO^^cm"^, a nd a power-law spect rum with r=1.73lQ 28. '^^^ following year a 
targeted radio search by [Roberts et al.l (120021 ) discovered J2021; these authors reanalyzed 
the ASCA data and concluded that a thermal component (kT ~0.1 keV) in addition to 
the absorbed power law improved the fit to the data, giving an = 7.Qt.t'l x lO^^cm"^ 



Hessels et al.l ( 120041 ) were able to resolve the inner equatorial PWN and jets surrounding 
J2021 with a 19 ks Chandra ACIS-S observation, noting a toroidal morphology and possi- 
ble double structure; this observation had a sub-framed field of view and could not probe 
the faint larger-scale structure. These authors report an Ah of 7.8]*]} 4 x lO^^cm"^, a power 
law index of r=1.73lQ2 ^^e extended emission. Thermal emission from the pulsar at 
kToo=0.15+0.02 keV was noted, along with a possible pulse fraction of 65% from the pulsar, 
presumably mostly due to thermal emission. We describe here a deeper Chandra exposure, 
which better resolves the innermost structure and maps the outer portions of the nebula. 
The improved constraints on spectral and morphological parameters of the PWN afforded 
by this deeper observation can help address the question of the pulsar's distance and will be 
helpful in interpreting future high energy observations. 
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2. Observations and Data Analysis 

We observed J2021 with CXO in the ACIS-S configuration (chips 2-3 and 5-8), with the 
pulsar positioned near the standard aim-point on the S3 chip during two epochs: December 
25-26 2006 with a total live time of 33.8 ks (obsID 8502 = 'obsl') and December 29-30 2006 
for 59.4 ks (obslD 7603 = 'obs2'). All figures shown below are merged images of the two 
observations. The CCDs were operated in VFAINT mode to improve rejection of particle 
backgrounds. Neither observation exhibited significant fiaring, so we include all the data in 
the analysis, for a total live time of 93.2 ks. These data are also compared with an archival 
19.0 ks observation on chips 6-8 taken on February 12 2003 (obslD 3901 = 'obsO'). For that 
exposure, the chips were ^ yindowed to 1/4 fra me, covering only the bright central region of 



the PWN, as described in iHessels et all ( l2004j ). 



The data were analyzed using CIAO 3.4.0 and CALDB 3.4.0. We started with level=l 
event files, removed pixel randomization, applied time-dependent gain and CTI corrections, 
and implemented a background cleaning algorithm for VFAINT mode which uses the outer 
16 pixels of the 5x5 event island to improve discrimination between good events and likely 
cosmic rays. HEASARC's WebPIMMS tool estimates a small pileup fraction of 1.6% at the 
pulsar location. Since the pileup is low, we improve the spatial resol ution of the ACIS image 



by applying an algorithm to correct the position of split pixel events iMori et al.l (120011 ): this 
decreases the on-axis PSF width in our data set by ~ 13%. 

In order to determine the background count rate, we located a region of the S3 chip 
relatively free of diffuse emission, cut out point sources, and measured a count rate of 0.59 
counts s~^ in the 0.5-7 keV band. Since much of the S3 chip is covered by diffuse nebulae, 
we split the background region into two parts: the southern corner of the chip, as well 
as a portion of the chip northwest of the pulsar, exterior to the extended emission. The 
nominal background count rate for the S3 chip in this band is 0.32 counts s~^. Although no 
large background flares were seen, solar activity was high during the observations which may 
contribute to the increased background. It is also likely that low surface-brightness diffuse 
emission from either extended PWN structure or Galactic background makes a non-negligible 
contribution. 



2.1. X-ray Spatial Analysis 

In order to bring out the large scale diffuse emission, we excised point sources (retaining 
the pulsar and a few other sources of interest), exposure-corrected the image to minimize 
chip gaps, and smoothed with a 6" FWHM Gaussian. Figure 1 reveals considerable emission 
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Fig. 1. — ACIS-S 1-7 keV image (exposure corrected, 6" Gaussian smoothing). Point sources 
have been removed except for the AGN candidate, a bright USNO star, and four sources 
coincident with a diffuse IR nebula. Extraction regions for the outer nebula (and dotted 
exclusion region) are also shown. 



surrounding the pulsar located at R.A. = 20:21:05.43, decl. = +36:51:04.63. A thin arc of 
emission extends I'.l east from the pulsar to the edge of the S2 chip. To the west, diffuse 
emission covers much of the S3 chip and extends onto S4. The southern boundary of this 
emission is very pronounced, with a sharp arc bracketing the pulsar to the east and south. 

The pulsar is surrounded by an axisymmetric inner nebula extending 20" x 10" with 
two ridges of emission flanking the pulsar point source. Southwest of the pulsar is a very 
bright inner jet 4" long, while to the northeast a similar, albeit dimmer, feature forms an 
inner counter jet. On larger scales the jet can be followed in a fainter extension from 4" 
to 10", and intensifies to a bright, narrow outer jet running 12" to 30" from the pulsar. 
This outer jet ends in a distinct knot of emission. On the opposite ('counter-jet') side the 
extended jet is not obvious. However, there is a hard, unresolved source 23" from the pulsar 
roughly along the counter-jet axis with no obvious stellar counterpart; this is plausibly the 
counterpart of the jet knot to the south. The emission from the outer jets may be variable. 



For example iHessels et al.l (l2004j ) did not conclusively detect the outer jet. In obsO, the 



structure is not visually apparent with 9.0±4.7 x 10 ^ counts s ^ from this region, while in 
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the new data the outer jet produces a significant detection with 1.3±0.2 x 10~^ counts s~^ 
above a surrounding bac kground annu l us. T here is also no detection of the jet terminal 
knot in obsO. In contrast. iHessels et al.l (120041 ) describe a diffuse structure on the counter-jet 
side. We remeasure this region and find 1.9±0.5 x 10~^ counts s~^ in obsO; the same region 
in the new data contributes 1.0±0.2 x 10~^ counts s~^ above the local background. While 
the excess is formally significant, there is no obvious coherent structure extending to the 
northern knot. Of course, variation in extended jet structure is not unexpec ted: the Vela 



pulsa r outer jet displays dramatic variability on timescales of days to weeks (jPavlov et al. 
2003h . 
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Fig. 2. — 1-7 keV images. Left: A soft stretch shows the double 'ridge' structure of the 
equatorial PWN, along with the extraction aperture of the outer jet. Right: A deeper 
stretch and 1" Gaussian smoothing bring out the diffuse emission and jets. The extraction 
apertures for the pulsar, inner jets and equatorial PWN (with dashed exclusion region) are 
shown against the PWN. 



Turning to the central, brightest portion of the PWN, we see two ridges of X-ray emission 
perpendicular to the jet symmetry axis, o ne crossing ju s t abov e the pulsar position and one 
crossing the bright inner jet. As noted in iHessels et al.l (120041 ). this suggests a double torus 
structure, like that of Vela, but seen more nearly edge-on. We further explore the PWN 
structure by fitting torus models in §2.2 below. There is a tradition of naming PWNe after 
animals, dating back to the fancied resemblance of the remnant of SN1054 in a drawing from 
visual observations by William Parsons, Earl of Rosse (1844) to, variously, a crab or a crab's 
claw. For the J2021 PWN, we have a double ridge extending to each side of the pulsar and a 
single bright narrow jet. We thus christen this the 'Dragonfly Nebula' with the double ridge 
forming the paired wings and the outer jet forming the tail. 
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As usual in any deep Galactic exposure, several dozen unresolved sources, many co- 
incident with field stars, are detected in our image. We mention here only the sources of 
particular interest. One of the brightest sources in the field, CXOU J202117. 4+364723. 7, 
has a hard, absorbed spectrum and no optical, IR, or radio counterpart. We consider this a 
likely AGN viewed through the plane. Roughly 4^7 southwest of the pulsar, a strong X-ray 
source is seen coincident with the bright field star USNO-Bl.O 1268-00448692. 

Examination of MSX 8/im data shows a cavity coincident with the X-ray nebula. In 
particular, the southern edge of the nebula coincides with a fairly sharp step in the 8/im flux 
(Figure 3). This can be interpreted as a dust deficit in the X-ray nebula interior, caused 
by evaporation of embedded dust grains. The trail also seems to correlate with a faint 
1.4GHz continuum in the DRAO GPS maps, but more sensitive imaging would be needed 
to confirm the reality of this structure. At the west end of this cavity the MSX 8yum data 
shows a peculiar double-lobed diffuse (1^7 x 2') object, centered at 20:20:36.1 +36:50:37 
(Figure 3). The center of this objec t is within 40 ^ ^ of t he Siaghz = 25.8 ± 4.1 mJy source 
NVSS J202036+365123. Intriguinglv fPavlov et all J2OO7I ) noted a double-lobed mid-IR/radio 
nebula of similar size near PSR B1823— 13 and its PWN, although the pulsar proper motion 
argues against association in this case. Four of the faint X-ray sources are coincident with the 
IR nebula. Three of the X-ray sources lie within the lobes, while the fourth is in the bisecting 
dark lane. All appear absorbed with virtually no counts below 1 keV; the association with 
the IR complex, if any, is not obvious. 




Fig. 3. — Left: ACIS image with 3" Gaussian smoothing showing the PWN, jets, and outer 
nebula trailing off to the west. Right: MSX A-band image (8.3 fim) overlaid with heavily 
smoothed contours from the X-ray data - a 8.3 fim deficit follows the outer PWN; faint X-ray 
point sources are found in the bright double-lobed IR source. 
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2.2. Spatial Modeling of the Inner Nebula 



In the introduction we noted that many PWNe display axisymmetric torus and jet 
structures formed from a relativistic wind with significant latitudinal variation . MHD sim- 



ulatio ns of such winds (IKomissarov fc LyubarskyI l2004l : iDel Zanna et al.l l2006l : iBucciantini 
20071 ) show that Doppler boosted arcs and jets can appear in synchrotron emissivity maps, 
often associated with lines-of-sight tangential to discontinuities in the fiow. For the brightest 
PWNe (i.e. Crab and Vela) the CXO images are sufficiently detailed for comparison with 
these models. For more typical sources, only the overall scale, orientation and Doppler en- 



hancement can be measured. T 



of pulsar phvs ics (INg fc Rornani 



lese parameters can be useful for probing several aspects 



20081 ). Even when the PWN structures are less than visu- 



ally stunning, iNg &: Romani showed that it is possible to extract robust values for 



the basic angles and scales by fitting a simplified Doppler boosted torus model to the CXO 
images. Following their methodology, we fitted the 'Dragonfiy Nebula' with a double torus 
model similar to that of the Vela PWN, including a point source and a uniform background. 
The fitting minimizes residuals using a Poisson-based likelihood function. We estimated the 
statistical errors by re-fitting Monte Carlo simulations of the best-fit model. The systematic 
errors are es timated by in t urn r e-fitting the data with the jets and point source regions 
removed; see iNg fc Romanil (120081 ) for details. 



For comparison, we also tried fitting with a single torus model; in this case, a large torus 
thickness, or 'blur', is required for a viable fit. This single torus model does a substantially 
poorer job of matching the inner structure near the pulsar. For example, if we exclude the 
pulsar and jets and collapse the counts along the symmetry axis there is a clear dip between 
the tori in the nebula brightness profile, 0.5-2'.'5 SW of the pulsar. In this region the single 
torus model departs from the data at the Ga level. However the large 'blur' in the single torus 
model allows it to fit the general excess of diffuse emission skirting the PWN better than 
the sharper two-torus model (see Figure 4). This allows the fft statistic to be of comparable 
quality. Thus with a free fft, the double torus model is preferred, but not deffnitively so. 
If an extra 'halo' component were added to the model, the preference for the double torus 
would be stronger, but this is not really required by the present data quality. 

Figure 4 shows the data in 0.5-8 keV band compared to the best-fft models, with the 
parameters listed in Table 1: the axis position angle \E', inclination angle torus radius 
i?, blur 5 for the torus proffle, post-shock ffow velocity /5, and separation d. The best-fft 
orientation parameters are consistent in both models, though R changes slightly due to 
the different blur. The new results also show good agreement with the previous ffts by 
Hessels et al.l ( l2004l ). but provide more precise measurements. The relatively small errors 
indicate that the ffts are robust and not sensitive to the presence of the jets. In summary. 
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the numerical fits confirm our visual impression that the double structure better matches 
the data. However, unless we allow extra model components to absorb the extended halo 
counts or fix the blur 6 at small values for both models, the fit statistic for the single torus 
is not dramatically worse. Note that the Poisson statistics employed here preclude a direct 
comparison of the two models with the F-test. 




Fig. 4. — 0.5-8 keV image compared to the best-fit double tori and single torus models. The 
dotted ellipse shows the region excised to estimate systematic errors. 



2.3. Spectral Analysis 

Wc used CI AO's Sherpa environment (version 3.4) and XSPEC (version 12.3.1) to fit 
the spectra of the pulsar and nebula. For consistency, we bin all spectra to a signal-to-noise 
ratio of three; as a result the jet was fit only with the new data sets, as the old observation 
(obsO) has insufficient counts. The extended emission lies outside the field of view of obsO, 
and therefore must also be fit with the new data sets alone. Careful background subtraction 
is important for robust spectral fitting; for the inner jets we subtract the surrounding inner 
nebula as background, for the star and AGN we utilized surrounding annuli, and for all 
other regions we used a large source-free background region. This background region is on 
the S3 chip for all spectral fits except for the eastern wisp, which uses a background region 
on the S2 chip. Fit errors are projected multi-dimensional values at the 90% confidence level, 
except for fluxes. Multi-dimensional flux error estimates arc often quite large due to spectral 
parameter uncertainties, so we list single parameter errors at the 90% confidence level. 

Wc determined the interstellar absorption for the pulsar/PWN complex by simultane- 
ously fitting the inner and outer components of the nebula with a global A^Hj but different 
power-laws; this yields a hydrogen column density of G.T^oiyXlO^^ cm~^, consistent (albeit 
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with much smaller er r ors) w ith the lRoberts. Romani. &: Kawai I (120011 ) v alue of 5.0±2. 5 x 10 



21 



Roberts et all (120021 ) value of A^^h = 7.6+^;^ x lO^^ cm'^ and the iHessels et all (120041 ) 



cm" 

value of 7.8l}'4xl0^^ cm~^. The value of 6.7 x 10^^ cm~^, is adopted in all final fits to com- 
ponents of the PSR/PWN complex. Although the soft thermal spectrum of the pulsar itself 
provides many photons, the effective temperature is highly covariant with N-h_; accordingly 
fits to the relatively simple power laws of the extended non-thermal emission provide the 
best absorption constraints. The two brightest point sources have absorptions bracketing 
the pulsar value. The field star Bl.O 1268-00448692 gave a low A^'h of 0.2 xlO^^cm^^. A 
power-law fit to our candidate AGN gives A^h = 2.2 xlO^^ cm~^, ~ 3x larger than the 
PSR/PWN value. This supports its identification as an extragalactic source, although some 
of the absorption may be intrinsic. 



2.3.1. Pulsar Spectrum 

To measure the spectrum of the point source with minimal nebular contamination, we 
define a source region of radius 1.5 pixels (0'.'74). Pulsar spectral fitting is complicated by 
the fact that the inner nebula surface brightness scales inversely with radius; extrapolating 
the radial profile indicates that the nebular surface brightness is ~ 2.5 x greater in the 
source region than in an annulus extending 2.5 to 6 pixels radially from the pulsar. This 
precludes simple annular background subtraction. Instead we determine the nebular power 
law in this surrounding annulus and scale up the amplitude as predicted from the brightness 
profile, giving an unabsorbed 0.5-8 keV flux of 4.1xl0~^^erg/cm^/s and r=1.20. We hold 
this amplitude and F fixed as a contribution to the background in the pulsar aperture, while 
fitting the pulsar point source spectral model. Both the pulsar and annulus regions utilized 
the large source-free background region on the S3 chip. 

For point source fluxing, we modeled the necessary aperture correction by simulating 
10 monochromatic PSFs from 0.5 to 9.5 keV with the Chandra Ray- Trace program, ChaRT. 
The enclosed energy fraction as a function of radius was then calculated to correct the 
ARFs used in the spectral fit. We fit four different models: a blackbody (BE), a power 
law plus blackbody (PL+BB), a magnetized neutron star hydrogen atrnosphe re (NSA), and 



a power law plus neutron star atmosphere (PL+NSA) (jZavlin et al.l Il996l ). all with the 
extra fixed power law component from the nebular background. As a check, we also fit 
allowing A'h to vary. The fitted model parameters were consistent with those found using 
the A'h fixed by the nebula fits, albeit with increased errors. For multi-component models, 
covariance between components will induce larger effective flux errors - e.g. for the power 
law plus atmosphere model the power law component flux errors increase from ~22% to 
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~26% allowing for component covariance (spectral parameters fixed). The errors grow to 
~100% if all spectral parameters are free. 

Adding a second thermal component to the simple blackbody model improves the fit, 
but since the pulsar has significant counts above 3keV, a high temperature T2 ~ 1.5 keV 
and small area A2 ~ 10~^yli would be required. A similar result is achieved for a two 
temperature neutron star atmosphere fit. The small area of the secondary component does 
not seem compatible with a polar cap for a Vela-type pulsar. The power law plus blackbody 
and power law plus neutron star atmosphere models also yield much better fits than the 
simple absorbed blackbody or NSA models. This point source non-thermal PL component 
may represent either magnetospheric emission or unresolved PWN structure. We favor the 
former, since the component contributes more flux than the nebular extrapolation. Although 
the nebular extrapolation is uncertain, varying over a plausible range only induced ~ 10% 
changes in the point source PL flux. 



Our PL-I-BB model returns a kToo of 0.16 ±0.02 keV which matches well to the Hessels et al 



(120041 ) fit of 0.15 ±0.02 keV. This value is also 



comparable to that of the similarly-aged Vela 



pulsar (r, = 11 kyr, kT^ =0.128±0.003 keV, 



kyr, kToo=0.17±0.02 keV, iRomani et al.l l2005l ) . The thermal fits allow us to calculate an 



Pavlov et all 120011 ) and B1706-44 (r^ = 17 



effective radius as a function of distance, and hence the emitting area. For the moment we 
set the scale with a nominal distance of lOkpc (so d = lOciiokpc), although as we see below 
smaller distances are preferred. Our fit to a power law plus blackbody with Nn fixed gives 
R = T.O^f^dio km, representing emission from 29t^2% of the stellar surface for a Roo = 13.1 
km star. The minimum emitting area of S.SxlO^^dfQ cm^ is far greater than that of the 
traditional polar cap Ape = 2tt'^R^/cP ~ 6.3 x 10^ cm^ for any reasonable distance. 

Blackbody deviations for NSA models in principle allow one to determine both radius 
and surface redshift, although in CCD quality data these are highly degenerate. Conse- 
quently, in fitting model atmospheres we fix the surface radius at Rs = 10 km, corresponding 
to Roo = Rs{l - 2GM/RsC^y^/^ = 13.1km. The neutron star mass is fixed at M = 1.4 M©, 
and we use models with a surface magnetic field of lO^^G. With the mass and radius fixed, 
the pulsar distance is determined directly from the model normalization. With A^h fixed 
our NSA±PL fit yields a distance of 2.1;J^^ q kpc, where the distance error includes multi- 
parameter uncertainty. As usual the NSA model's Wien excess requires a lower Teg and 
a smaller dis t ance than that found from blackbody model fits. For a similar NSA±PL fit 



Pavlov et al.l (120011 ) report a Too of 0.059 keV from Chandra ACIS-S observations of the Vela 
pulsar. These authors infer a distance of 220±20spc from their neutron star normalization, 
which is within 25% of the accepted ~ 300 pc distance. Our J2021 neutron star atmosphere 
normalization distance estimate of 2.1 kpc therefore lends some credence to the assertion 



that the pulsar is significantly closer than 12 kpc. 



2.3.2. Extended Source Spectra 

As a check, we fit the spectra for the inner and outer nebulae independently, allowing 
A'h to vary; we find A^h = 6.7^} o x 10^-^ cm~^ and A'h = 6.6;^]]; 'J x 10^^ cm~^ respectively, 
fully consistent with the global absorption fit above. Using the fixed global A^h allows us to 
extract useful spectral index estimates for various components of the PWN complex. The 
large arc east of the pulsar may or may not be associated with J2021. Accordingly we also 
fit this spectrum allowing A^h to vary, though with the larger errors the result is not very 
constraining. 

As one moves progressively further from the pulsar the extended emission exhibits sig- 
nificant softening, as expected from synchrotron burn-off in the outer PWN. The inner jets 
have the hardest spectra (r=l.l), while the outer jet region is softer with r=1.7. The inner 
nebula or 'equatorial PWN,' has P = 1.5. The diffuse trail close to the pulsar («i30"-60") 
has a power- law index of 1.7, while in the farther reaches of the outer nebula we observe an 
index of 1.9. 



2.3.3. USNO Star Spectrum 



We estimated a photometric parallax for USNO-Bl.O 1268-00448692 using data from 
the USNO-B (~B,R,I bands) and 2MASS catalogs (J,H,Ks band s) . We converted the survey 
magnitudes to standard colors an d fit the colors to ste llar values ( iKenyon &: Hartmanrull995l ) 
subject to interstellar reddening (ISchlegel et al.lll998l ). The best global fit to main sequence 
colors implied class M2 V, with Ay ~ 0.45 at a distance of ~ 60pc. The POSS-II 'P band 
point has the largest departure from the model colors; if it is excluded the best class shifts 
to K6 V with Ay ~ 1.7 at ~ lOOpc. Few post-MS stars of this color emit coronal X-rays. 
The best match (excluding I) was a rather implausible G5 III with Ay ~ 2.7 at 1.5 kpc. 

The star sh ows modest va r iabilit y within and between epochs. We therefore follow 
the procedure of iGetman et al.l (120051 ) and adopt a two component Mekal plasma model; 
we fit the temperatures and amplitudes separately for each observation, while fitting for 
a common A^h = 0.2 x lO^^cm"^. As expected the star temperatures are higher during 
the X-ray bright epochs. The fit Ah can be compared with the optical extinction, where 
A^H = Ayl.Q X lO^^cm"^ predicts Ay = 0.4:~^q^. The best global interpretation thus appears 
to be a relatively mundane nearby M2 V field star; the estimated A^h and Ay are in reasonable 
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agreement, albeit somewhat higher than one would expect for a 60 pc distance. However this 
extinction adds no significant contribution to the full column to J2021. 



Discussion 



3.1. Distance 



As noted in the introduction, the dispersion meas ure of 369cm~^p c (A^p,= 1.13xl0^^cm~^) 
for J2021 corresponds to a distance of 12kpc in the ICordes &: Lazid (|2002| ) model. Such a 
large distance is problematic on several grounds, so it is important to check this estimate 
against other observables. Our best fit hydrogen column density is A^H=6.7xl0^^cm~^, while 
the total Gal actic A^h in the direction o f J2021 is estimated from HI m aps to be 1.2xl0^^cm~ ^ 
according to lDickey &: LockmanI (jl990l ) or 9.7xl0^^cm~^ according to lKalberla et al.l (120051 ). 
Thus we already detect somewhat less than the expected column, suggesting that J2021 does 
not reside in the far reaches of the Galaxy. Moreover, with Ah = 2.2 ± 0.9 x lO^^cm"^ mea- 
sured for our AGN candidate an even higher column is implicated (although, again, some of 
this absorption may be intrinsic). On the balance, these data suggest that the pulsar/PWN 
absorption, while large, is ~ 1/3 — 2/3 the full Galactic column density in this direction, 
placing the source well inside the Milky way disk. 



According to lCordes fc Lazid (120021 ). DM distance estimates are seldom off by > 50%, 
unless the line of sight intersects anomalously dense electron clouds. In any case, at such large 
DM it is difficult for any one HII region to substantially increase the value. Nevertheless, 
we note that our observed n eutral to ioni z ed rat io An/A'e ~ 6 is somewhat below the 
canonical factor of 10, while iHessels et al.l (120041 ) noted that J202rs radio pulse profile 
showed anomalous scattering with a scattering measure SM ~ 100 x that predicted by the 
NE2001 model. These provide evidence for high ionization along the line of sight. However, 
none of the X-ray structures that we see provide an obvious source of such ionized gas: the 
'Arc' passes north of the pulsar and the extended nebula is too faint to represent a large 
emission measure or column density. 

One problem with the large 12kpc distance would be the high 7-ray efficiency. For an 
assumed Isr beaming, the observed 3EG point source 7-ray fiux gives a 0.1-5GeV efficiency 
of 7]^ = L^/E ~ 0.15diQ. We should compare this with the similar 7-ray pulsars Vela 
(d^ 300pc, E = 6.9 xlO^^ ergss'^) and PSR B1706-44 (d~3 kpc, E = 3.4 xlO^*^ ergss^^) 
which have Isr efficiencies of rj = 6.6 x 10~^ and 1.8 x 10~^, respectively, for the same energy 
range and assumed spectral index. If the efficiency matches that of PSR B1706— 44 with 
nearly identical spin properties, we infer a distance of ~ 4kpc. At the very low efficiency 
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inferred for the Vela pulsar, the source would need to be at a rather implausible 800 pc. 
We can also use trends in the X-ray luminosities of ot her pulsars and PWNe to obtain 



rough distance estimates. For 6 other Vela-like pulsars iKargaltsev et al.l (120071 ) find an 
efficiency for conversion of spindown power ranging from ?7x ~ 3 x 10^^ to 3 x 10^'^ in the 
0.5-8 keV band. Some of this spread likely stems from variation in the ambient medium as 
well as differing pulsar properties. If we adopt only the unabsorbed flux of the equatorial 
PWN and the jets, we see that the observed 5 x 10~^^ergscm~^ s~^ implies a distance range 
of 1.3-4.1 kpc. Inclusion of emission from the extended outer nebula does not seem consistent 
with the flux estimates for other PWNe, but would shrink the implied distance range to 0.7- 
2.3 kp c. The non-therraal lura inosity of the point source should also scale with spin-down 
power; IKargaltsev et al.l (120071 ) note that these quantities correlate surprisingly well, with 



^PWN ~ 5Lp°"*^. For J2021, the power law flux is relatively poorly determined, but at 
only 0.5 x 10"^'^ergs cm~^ s~^ even the equatorial nebula represents ~ 8x this component, 
making it easiest to interpret the equatorial nebula as the full PWN flux and supporting the 
1.3-4.1 kpc distance. We caution that this estimate should not be taken too seriously given 
the large scatter in efficiency ratios. 

A final spectral estimate of the distance comes from the neutron star thermal flux, which 
if interpreted as magnetic H atmosphere emission from the full surface of a canonical neutron 
star, implies d ~ 2.1i?iokpc, for a 10 km local neutron star radius. If the surface spectrum is 
a blackbody plus power law, then larger distances (4-8 kpc) are nominally allowed. However 
the required ~ 2MK temperature would be anomalously high for full surface emission of a 
Vela-aged pulsar. A polar cap could be this hot, but an implausibly small ~ 0.8 kpc distance 
would be required to match the flux to the canonical polar cap area. In sum, though, the 
spectrum and flux estimates for the X-ray emission from this pulsar and it's PWN, along with 
the comparison of the inferred 7-ray efficiency with other Vela-like pulsars implies a distance 
of 3-4 kpc, substantially smaller than the DM-inferred distance. We adopt an estimate of 
d = Adi kpc , with ^4 1 in what follows. The source of the extra dispersion and scattering 
along the line-of-sight to J2021 is yet to be identified. 



3.2. Birthsite and Environment 

Most pulsars show bow shock PWNe when traveling at supersonic speeds through the 
interstellar medium. However, with a characteristic age Tc = 1.7 x lO^yr, J2021 should still 
reside in its parent supernova remnant, where in the absence of a central pulsar or PWN the 
reverse shock reaches the center on a timescale of 

tsed ~ 7{M,,/10MQf\EsN/lO''ergs)-'/\njsM/lcm-Y'/'kjT (1) 
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( iReynolds fc Chevalierl Il984l ). In the presence of a young pulsar, the reverse shock will 
collide with the expanding PWN in a time somewhat less than this. If the remnant is just 
entering the Sedov phase the reverse shock should be reaching the remnant center, where it 
would crush the PWN. For typical high pulsar velocities the system can be well off-center in 
the parent SNR; the crushed P WN should then trail back toward the sup ernova site, with 



the pulsar at the leading edge (Ivan der Swaluw. Downes fc Keeganll2004l ). Instead, we see 



symmetric equatorial structure, with jets extending to ~ 30" and no prominent bow shock 
on smaller scales. This suggests a low velocity pulsar in a relatively quiescent surrounding 
medium. It is possible that the reverse shock has not yet reached the PWN or, alternatively, 
that the reverse shock reverberations settled down some time ago, leaving a relatively uniform 
high-pressure SNR interior as the ambient environment. Given the characteristic age of the 
pulsar (17 kyr), and the timescale for the SNR reverse shock to reach the PWN (~ 5 kyr), 
we favor the latter interpretation. 

If a lO^^E'sierg supernova has produced a Sedov-phase remnant of age 10^t4y in an 
external medium of density riisu cm~'^, then we expect a shell of radius 

^SNR ~ l2'{E^JnisMf'hTld^. (2) 

While there is no cataloged remnant and no obvious X-ray or radio shell in the vicinity, we 
could plausibly associate the thin arc with a partial SNR shell. However, with a radius of 
curvature ~ 7' centered ~ 8' northeast of the pulsar, it seems unlikely that this represents 
the limb of the parent SNR. Of course, the scale would be acceptable for the large ~ 10 kpc 
distance, but then the pulsar would have left the shell at 2300(iio/t4 kms~^; the PWN should 
then show a bow shock with an unresolved stand-off. It is possible that the arc is a filament 
on the face of a larger remnant, with the pulsar close to the center. A search for a large 
~ 30' low surface brightness shell could be productive. 

Assuming for the moment that the PWN does live in the heated interior of a Sedov- 
phase SNR, we expect an ambien t pressure Psnr pR^ = xp^sm-^sn^'^^^ ■ ^^^^ X = 



0.047 (Ivan der Swaluw et all 120001 ) so that 



SNR 



lxlO-'El(\f^,,tf/'gcm-'s-' (3) 



The PWN equatorial torus represents the termination shock where the pulsar wind momen- 
tum flux balances this pressure. Two factors complicate our estimation of the termination 
shock radius for J2021. First the wind is expected to be equatorially confined so that the 
momentum balance condition is 

PsNR ~ E^/iAnRl^sc) (4) 
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with ^ > 1. For example, for a wind momentum scaling with latitude 6 as sin^ (iBogovalov fc Khangoulyan 



20021 ) one expects ^ = 3/2. Second, while in some PWNe such as the Crab, we observe the 
'sub-luminous zone' of unshocked pulsar wind and hence directly resolve Rws^ here we only 
see the overall extent of the equatorial flow, which is generally 2-3 x larger (e.g. Ng & Ro- 
mani 2008). Accordingly, using the pressure (3) we estimate the overall angular scale of the 
equatorial torus as 

^ (2 - ?,x)V!QE-^'\-l'^hTe''d-,\ (5) 

Taking E^i = rifSM = 1, ~ 2 and t4 = 1.7 (the spindown age) we get an overall scale 
of 9eq ~ (2 — 3x)3.1"d^^ . If we identify the double ridge of the 'wings' of the dragonfly with 
this overall torus, we have a measured size of ~ 10", which from (5) implies a distance of 
2.5-3.7kpc. This is in good accord with our previous estimate of 3-4 kpc. 

More heuristically, we can simply scale to the observed angular radius of the overall 
toroidal flow for Vela (~ 30") and PSR B1706-44 (~ 15"). Matching the ~ 10" of the 
dragonfly wings suggests ~ 1 or 4.5 kpc (for Vela and B1706, respectively). Again, while not 
precise, the geometrical estimate of the wind shock scale suggests a distance < 4 kpc. 



3.3. Outer Nebula and Jets 

In Figure 4, we see that the diffuse emission surrounding the equatorial PWN is brightest 
in a ~ 25" radius halo around the pulsar, continuing to the 'outer nebula' which forms a 
~ 50" radius limb-brightened structure surrounding the pulsar and trailing off to the west, 
with some curvature. If we interpret this outer nebula as a bow shock, then with a PWN 
speed of lOugkms"^ {v = vq x lO^cms"^) relative to an ambient medium of number density 
n, we have 

Obs ^ [E/iAncpvlsnd')]'^' ^ AQ"n~'" / {v,d^) (6) 

and the observed O^s ~ 50" standoff distance implies a rather low speed of vpsR ~ Sn^^^'^d^^ 
kms~^. This nebula seems to extend ~ 4' to the west. If this traces the actual motion since 
pulsar birth we require a fairly typical pulsar speed of f = A'd/Tc = 270(^4 km s~^ assuming 
the true age is the spindown age. To preserve the large nebula standoff, the ambient medium 
should then have a density n < lO^^cm^^, which is plausible for the hot interior of a SNR. 
Some support for this interpretation comes from the softening of the X-ray spectrum in the 
western half of the nebula, which could be attributed to synchrotron cooling of relic particles 
at the pulsar birth site. Alternatively, this outer nebula could be caused by PWNe electrons 
filling a pre-existing structure in the SNR interior, e.g. density variations in a progenitor 
wind, so it might not trace the pulsar motion. 



- 16 - 



Turning to the jet components, it is natural to attribute the flux ratio between the jet 
and counter jet to Doppler boosting. For a continuous co-hnear inner jet /counter-jet we 
expect Ij/Icj = [(1 — /5cosC)/(l + /5cosC)]~^~\ where ( is the inchnation of the jet axis to 
the line of sight. The spectral fits give F ?a 1.1 and a flux ratio of 1.7]^q5. However, the 
spatial flt to the the equatorial structure places it nearly edge-on {( ^ 85°), which gives a 
maximum flux ratio of ~ 1.4 for large /3. Of course even a small jet misalignment relaxes 
this stricture; for ( ^ 80°, we can accommodate a flux ratio of 1.7 with (3 = 0.7. Still, with 
the jets nearly in the plane of the sky large f3 are required for signiflcant Doppler boosting. 

A similar argument applies to the outer jets. For the jet an aperture that avoids the 
bright knot at the end provides 73 ± 16 counts, while the same size aperture for the counter- 
jet has 28 ± 14 counts (not a signiflcant detection). This leads to a fairly large flux ratio of 
2.6^09. Of course this cannot be accommodated for an inclination as large as 85°. However, 
visually the jet and counter-jet do appear to be mis- aligned on the plane of the sky, so a 
larger counter-jet inclination, allowing larger Doppler de-boosting, does seem plausible. For 
example, for the observed jet spectrum (F ^ 1.7) any angle < 75° can be accommodated for 
(3 < 0.7. 

In view of the possible jet variability noted above, none of these flux ratio constraints 
should be taken too seriously. However, if we do attribute the faintness of the outer counter- 
jet to Doppler boosting, then the presence of a terminal knot for both the jet and counter-jet 
could be due to a termination shock, resulting in slower flow with more isotropic emission. 
The outer nebula is somewhat brighter in a ~ 25" region around the PWN and perhaps 
this marks the flow boundary where the jets terminate. This would be reasonable if they 
lie close to the plane of the sky. In any case the termination probably lies well within the 
envelope of the 'outer nebula;' with the jets only half the length of the outer nebula radius, 
an implausibly small inclination of 30° would be required to make their projected length as 
small as observed. 



3.4. Energetics and Spectral Results 



The observed spectral softening from the inner to the outer nebula (Figure 5) suggests 
aging of the synchrotron population on the flow time. At radio frequencies, however, the 
dominant population should be uncooled with a spectral index (Xr ^ 0.3. An exarn ination 
of the 1.4GHz DRAO Galactic Plane Survey continuum maps (lEnglish et al.l Il998l ) shows 
a 3±0.5mJy source at the position of J2021, not clearly resolved at the ~arcmin GPS 
resolution. Since the pulsar itself provides only 0.1 mJy at 1.4 GHz, we interpret this as 
PWN emission. 
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Fig. 5. — PWN synchrotron cooling. Here we show the radial softening of the PWN X-ray 
spectrum (see Table 2). For the dominant flux of the outer nebula we extrapolate to the 
observed radio flux with a Aq;=0.5 cooling break. 



For such a low radio flux, the cooling break must be at relatively high frequency. Ex- 
trapolation of the X-ray flux with a Aa = 0.5 break imphes Log[i/B(Hz)] = 14.9lQ;g. The 
synchrotron break frequency at an age of 10^t4yr in a field of SmG^iG is z/^ = 1.0 x 
10^°-B~^t4 ^ Hz. Thus we infer an outer nebula field of i? = 13^4^ ;uG for t/^ = 1.7. This can be 
compared with the field extrapolated from Bg = 3.2 x lO^^G at the pulsar surface (7?s=10km). 
If the field is dipolar within the magnetosphere, followed by cx 1/r in the wind exterior to 
the light cylinder {Rlc = c/lj) we expect By^g ~ ^BsR^/ {R\cRws) = 2206^ ^d4^^iJ,G in the 
wind termination shock at an angle of 59'^. If the field continues to decrease as oc 1/r in the 
post-shock flow, then we would expect a value ~ 20//G at the ~ 50" limb of the outer nebula. 
This is in reasonable accord with the photon- weighted value estimated from the (rather un- 
certain) spectral break. The corresponding magnetic pressure B'^/8n ~ 2 x 10~''^"^gcm~^s~^ 
is somewhat lower than equipartition for a SNR interior. More detailed radio measurements, 
resolving the PWN and determining the radio spectral index would substantially tighten 
these arguments. It may also be useful to mount deep optical/IR imaging to detect the 
PWN flux near the cooling break. 
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Conclusions 



We observed PSR J202H-3651 and its surrounding structured nebula and jets with 
the CXO ACIS. This 'Dragonfly Nebula' displays an axisymmetric morphology, with bright 
inner jets, an apparently double-ridged equatorial inner nebula, and a ~ 30" long polar 
jet. Surrounding the central nebula is a low surface-brightness outer nebula. This structure 
brackets the pulsar at a radius of ~ 50" and trails off to the west over 3-4'. The overall 
str ucture is highly rem iniscent of CX(9/ACIS image of the PWN surrounding PSR B1706- 
44 (IRomani et al.ll2005l ). and, as for this source, is plausibly caused by a low- velocity pulsar 
in a relatively quiescent medium. No clear evidence of a host SNR is seen. IR data does 
appear to show a dust deflcit associated with the trailed X-ray nebula. 

Although the photon statistics are limited, double structure for the equatorial PWN is 
suggested by both the initial 19 ks exposure and, independently, by our new deeper image. So 
far Vela is the only other pulsar showing a double ridge in its PWN. Whether this represents 
a physical doubling of the equatorial torus or projection (caustic) effects in an optically thin, 
Doppler boosted pulsar wind is still unclear. 

The new spectral data combined with the spatial scale of the PWN shock and its 
surrounding nebula do not support the ~ 12 kpc distance suggested by the large pulsar DM. 
The X-ray absorption, efficiency arguments, and comparison with other Vela-type PWNe, 
while individually inconclusive, all suggest a smaller d ~ 3 — 4 kpc. On the balance we prefer 
this distance scale, although we find no obvious source of ionized plasma that could explain 
an anomalously large DM, so this large value would remain unexplained. 

This new distance estimate will likely be helpful in interpreting new 7-ray observations. 
At this distance, the GeV efficiency needed to produce the EGRET source is comparable to 
that for PSR B1706— 44. Upcoming AGILE, and especially GLAST, observations should test 
this association. A high quality pulse profile will also be very useful, as the pulse shape is 
sensitive to the inclination angle C estimated here from the PW N fit. For example, in outer 
magnetosphere gap models (e.g. iRomani fc Yadigaroglul Il995l ). the large inclination angle 
of this pulsar predicts a double peaked 7-ray light curve with large phase separation. At 
higher TeV energies, the PWN may also be detectable, although the striking morphological 
similarity to the PWN of B1706— 44 inspires a note of c aution, since H.E.S.S. searches for 
B1706— 44's PWN find only upper limits of 0.01 x Crab (lAharonian et al.ll2005l ). However, 
it may be rel evant that Milagr o finds a source MGRO J2019-t-37 ~ 20 ± 25' west of PSR 
J2021-I-3651 (lAbdo et al.l 120071 ). Intriguingly, the offset from J2021 is along the direction 
of the 'outer nebula' trail. Since the TeV ICS emission is dominated by lower energy, and 
hence older, electrons than those that produce the X-ray synchrotron emission, such offsets 
are commonly seen for PWNe. A high quality lACT image from VERITAS or MAGIC 
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detecting the PWN and showing its arcminute-scale structure would be of particular interest. 
Conversely, strong upper limits on associated TeV emission would suggest that some common 
attribute makes the Compton component of the 'Dragonfly Nebula' and PSR B1706— 44's 
PWN anomalously faint. 
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double tori 



single torus 



cn 

Ri") 
blur (") 

d{") 

Point Source (cts) 
Torus (cts) 



50.1 ±0.4 ±0.6 
85.9 ±0.2 ± 1.0 
11.1 ±0.2 ± 1.0 
1.2 (fixed) 
0.84 ±0.01 ±0.001 
3.7 ±0.04 ±0.12 
1240 
3896 



50.0 ±1.1 ±0.11 
84.9l[5:^ ± 0.2 
9.3l°i±1.0 
2.9 ±0.04 ±0.2 
0.78 ±0.01 ±0.03 



1209 
4459 



Table 1: Best-fit torus parameters for the double tori and single torus models, with la 
statistical and systematic errors. 



Table 2: Spectral Fits to PSR J2021±3651 
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tat a distance of lOkpc 

* yields a distance of 2.1 kpc. 
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Table 3: Spectral Fits to Other Sources 
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AGN 


99 + 10 


1.7* 


0. 15 ± 0.07 


0.85 ± 0.11 


9.5/18 


1-2 


Star 


'-'•^^-0.20 


3.45lH?/16.ll^:r 


1.34 ±0.17 


1.43 ±0.19 


37.1/51 





Star 


n 20+°-^° 

U-^U_o.20 


3.8i;°:|/i5.8itf 


1.23 ±0.13 


1.31 ±0.14 


57.2/68 


1 


Star 


'-'■^'-'-0.20 


2.801^^/10.11^:^2* 


0.70 ± 0.08 


0.76 ± 0.08 


57.3/63 


2 



*held fixed 

■'■T(MK) for two Mekal thermal components 
^■0.5-8keV fluxes in units of 10"" erg cm s 



